Electrostatically defined few-electron double quantum dot in silicon 
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A few-electron double quantum dot was fabricated using metal-oxide-semiconductor(MOS)- 
compatible technology and low-temperature transport measurements were performed to study the 
energy spectrum of the device. The double dot structure is electrically tunable, enabling the inter- 
dot coupling to be adjusted over a wide range, as observed in the charge stability diagram. Resonant 
single-electron tunneling through ground and excited states of the double dot was clearly observed 
in bias spectroscopy measurements. 
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Electrostatically defined single and double quantum 
dot (DQD) systems in GaAs/AlGaAs heterostructures [TJ 
[2] are the current benchmark for the implementation of 
DiVincenzo's criteria using semiconductor qubits [3l[4l|5]. 
Although the nuclear spins inherently present in GaAs 
provide a fast decoherence mechanism, this drawback has 
been partly overcome recently [6 . Silicon has a natu- 
ral advantage in this respect since the only stable iso- 
tope with a nuclear spin is ^^Si. The 4.7% abundance 
of this isotope in "^*Si can be reduced by isotopic pu- 
rification, resulting in nearly nuclear-spin-free crystals. 
This should, in principle, increase the coherence time of 
electron-spin qubits in silicon [71 [8] . Initial demonstra- 
tions of Si-based DQD systems for spin qubits [10] 
have stimulated a number of recent studies of DQDs in 
both multi-gated silicon-on-insulator (SOI) [11] [12] and 
Si/SiGe [13j structures. 

In this letter, we report the fabrication of a few- 
electron DQD and its electrical measurement at milli- 
kelvin temperatures. The double dot is based upon a re- 
cently developed double-gated silicon quantum dot [M], 
which was also shown to operate effectively as a radio- 
frequency single electron transistor [T5^. Our approach 
provides a simple method of producing multi-gated sili- 
con quantum dots without the need for complement ary- 
MOS (CMOS) process technologies, such as polysilicon 
deposition and etching. The morphology of the double 
dot device is investigated using cross-sectional transmis- 
sion electron microscopy (XTEM) analysis. Transport 
spectroscopy demonstrates the ability to tune the double 
dot from the weakly- coupled to strongly-coupled regime. 
In the weakly-coupled regime, extracted capacitances of 
the system show good quantitative agreement with sim- 
ple modelling using FastCap [16]. 

The devices investigated in this work were fabricated 
on near-intrinsic silicon wafers (p > 10 kl]-cm at 300 K). 
After definition of n+ ohmic contacts by phosphorus 
diffusion through a masked sacrificial thermal oxide, a 
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FIG. 1: (a) SEM image of the Si MOS DQD. The three bar- 
rier gates and the top gate have widths ^30 nm and ^50 nm 
respectively. The Al barrier gates were plasma-oxidized to iso- 
late them from the top gate, (b) Schematic cross-section of 
the device (not to scale). Source and drain n^ contacts (red) 
were formed by phosphorus diffusion into the Si substrate 
(light blue). The top gate induces a 2DEG and the barrier 
gates create three potential barriers, forming two dots. The 
size of the dots is estimated to be 30 x 50 nm^. (c) Color- 
enhanced XTEM image of a similar device, (d) Enlarged 
XTEM image, showing sharp interfaces between the Si sub- 
strate, Si02 gate oxide, A\xOy and the Al top gate. 



200 nm field oxide was grown. In the active device re- 
gion (30x30 /im^), the field oxide was etched locally and 
replaced by an 8 nm-thick high-quality Si02 gate oxide, 
grown in an ultra-dry oxidation furnace at 800 °C in O2 
and dichloroethylene. Three Al barrier gates were then 
patterned by electron beam lithography (EBL), thermal 
evaporation and lift-off. The barrier gates were next 
passivated by plasma oxidation [Ml |17] , resulting in an 
electrically-insulating A^O^ layer surrounding the bar- 
rier gates. The Al top gate was defined in a second 
EBL step aligned to the lower gates with an accuracy 
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of ~20 nm. Finally, the devices were annealed at 400 °C 
for 15 mins in forming gas (95% N2/5% H2) to reduce 
the Si/Si02 interface trap density (Z^^t). Deep- level tran- 
sient spectroscopy of similarly-processed structures re- 
vealed Dit of order 5x10^^ cm~^eV~^ near the conduc- 
tion band edge [18]. 

Figures l(a,b) show a scanning electron microscope 
(SEM) image and a schematic cross-section of a dou- 
ble dot device. The top gate, which extends over the 
source and drain n+ contacts (not shown) and also the 
three barrier gates are used to form a two-dimensional 
electron gas (2DEG) accumulation layer under the thin 
Si02 layer. The barrier gates are used to locally deplete 
the 2DEG, forming three tunnel barriers that define two 
dots in series. The dots are geometrically defined by the 
distance between adjacent barrier gates (~30 nm), and 
by the top gate width (~50 nm). The outer barrier gates 
and top gate are used to control the electron occupancies 
electrostatically and the middle barrier gate is used to 
control the inter-dot coupling. 

Figure 1(c) shows an XTEM image along the top gate 
(i.e. perpendicular to the barrier gates). Apart from 
an increased (200 nm) top-gate width in order to aid 
XTEM sample preparation, this device is nominally iden- 
tical to the device used in electrical measurements. The 
XTEM image confirms the target 5 nm Al^^O^ layer thick- 
ness from the plasma oxidation process used (100 mTorr, 
50 W incident RF O2 plasma, 150 °C for 3 mins). Inte- 
restingly, at the interface between the top gate and the 
Si02, we find an additional Al^^O^y layer (~2 nm thick, 
see Fig. 1(d)) which could be due to the oxidation of the 
Al top gate via chemical interaction with the Si02 layer 
below. We note that the Al barrier gates, initially evap- 
orated to a thickness of 30 nm, show an Al core of only 
^20 nm in diameter after plasma oxidation, consistent 
with the formation of a ~5 nm A^O^y insulator. This 
A^O^ thickness is sufficient to allow differential biases 
of up to 4 V between the upper and lower gates with 
negligible leakage. 

Electrical (dc) transport measurements were per- 
formed in an Oxford Instruments Kelvinox KlOO dilution 
refrigerator at a base temperature of ~50 mK. A source- 
drain excitation voltage Vsd=^0 jiV at a modulation fre- 
quency of 13 Hz was used to monitor the differential con- 
ductance dl/dVsd- The source-drain dc current Isd was 
measured with a room-temperature current preampli- 
fier. Initially, the left (right) dot was characterized inde- 
pendently by setting the right (left) barrier-gate voltage 
Vbr{Vbl) equal to the top gate voltage Vtop- The mid- 
dle barrier gate voltage Vbm was fixed at yeM=0.818 V. 
Under these conditions Coulomb diamonds were recorded 
and the charging energy of the left (right) dot, was de- 
termined to be Ec '^5 meV (~2.5 meV) at Vtop =1-6 V. 
Therefore, the total capacitance of the left (right) dot was 
Cj:,ieftiright)=eyEc -30 aF (-60 aF) at Vbl=0.76 V 
{Vbr=0.76 V). To compare these experimentally ob- 
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FIG. 2: Differential conductance dl/dVsd as a function of the 
barrier gate voltages, Vbl and Vbr, for Vtop =1-6 V and 
zero source-drain bias. Tuning the middle gate voltage in the 
range Vbm=0.S14 V — 0.830 V, we observe a transition from 
two almost isolated dots (a) to the formation of a single large 
dot (b). 



tained results with modelled parameters, we used Fast- 
Cap which calculates the capacitances based on a finite 
element approach. Using the lithographic device dimen- 
sions as inputs, we obtained a total capacitance Cy: ~30 
aF for both dots, in good agreement with the experi- 
mental value for the left dot but at variance with that 
of the right dot by a factor of two. Such variations in 
capacitance from dot to dot could result from physical 
asymmetries in real devices, as evidenced by the XTEM 
image in Fig. 1(c), or from the presence of fixed charge 
in the gate oxide or at interfaces which can modify the 
effective gate potentials. 

We estimate the electron occupancy of a single dot 
using two methods. The first method uses Hall mea- 
surements of a similar MOSFET device from which the 
electron density is determined to be n = 3.5 x (Vtop — 
VTi/)10^^ cm~^ [19], where Vth is the threshold voltage. 
When operated as a simple MOSFET, our device showed 
Vth ~1.25 V. Hence, at Vtop=1-6 V we estimate the 
2DEG density of our device to be n ~ 1.2 x 10^^ cm~^, 
resulting in a dot occupancy of N~20 electrons for a 
30x50 nm dot size. As an alternative method, we es- 
timate electron occupancy by counting Coulomb oscilla- 
tions from Vth^ assuming no free electrons in the dots 
below threshold voltage [14]. This method derives a dot 
occupancy of N~15, in reasonable agreement with the 
previous method. Both approaches indicate that the de- 
vice operates in the few electron regime. 

Figure 2 shows the differential conductance dl/dVsd 
of the DQD as a function of the barrier-gate voltages, 
Vbl and Vbr^ for a fixed top-gate voltage Vtop =1-6 V 
and source drain voltage Vsd=^ V for two different mid- 
dle barrier-gate voltages Vbm- In Fig- 2(a), the rela- 
tively low middle barrier-gate voltage yeM=0.814 V and 
therefore high central barrier separates the two dots, re- 
sulting in the characteristic honeycomb-shaped charge 
stability diagram. By calculating the voltage ratios 
AVg^/AVBi? (AV^/AVbl), we can estimate the ratios 
of the mutual capacitance to the total dot capacitance 
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Cm/C^,ieft{right) ~0. 10(0.07), indicating that the double 
dot is in the weak couphng regime [2 . There, we observe 
the characteristic triple points resulting from the align- 
ment of the electrochemical potentials of the dots and the 
leads. In addition, current is observed along the sides of 
the hexagons, which can occur when the dots are strongly 
coupled to the leads and second-order co-tunneling pro- 
cesses occur [20] . Increasing the middle barrier gate volt- 
age to VeM =0.838 V, the mutual capacitance increases 
and dominates the system {Cm/C^^ieft(right) "^1)- This 
situation occurs when the middle barrier is reduced and a 
single (merged) large dot is formed, resulting in diagonal 
parallel Coulomb lines, as observed in Fig. 2(b). 

Figure 3(a) shows transport data through the DQD 
in the weak coupling regime Vbm =0.802 V with 
^TOP=1.4 V and Vsd=-1'0 mV. For \Vsd\ >0 the triple 
points evolve into so-called bias triangles, reflecting the 
occurrence of transport within the bias windows [2! • In a 
double dot system, two types of coupling can be distin- 
guished: capacitive coupling; and tunnel coupling. While 
capacitive coupling is a purely classical effect, tunnel cou- 
pling arises from the overlap of electron wave functions, 
classified by the fractional splitting ratio F = 2AVs/Vp^ 
where AVs is the splitting between the paired triangles 
and AVp is the diagonal separation between triangle pairs 
in Fig. 3(a) [2T1[^. Here, we find F -0.2, indicating that 
the interaction between the two dots is dominated by ca- 
pacitive coupling. The device may therefore be modelled 
using a capacitive approach. 

From the dimensions of the hexagon and triangles in 
Fig. 3(a) we obtain the key capacitances defining the 
system [2], namely: the total capacitances of the left 
and right dots, C^^ieftiright)'') the mutual capacitance be- 
tween the two dots. Cm] the relative capacitances be- 
tween each side barrier gate and its immediate neigh- 
boring dot, CBL{BR),ieft{right)] and the cross capacitance 
between each side barrier gate and the next neighbour- 
ing dot, Cy,BLixBR),rightiieft)' Thcsc rcsults agree weh 
with modelling performed using FastCap (see Table I). 
We note that by appropriate tuning of the barrier gate 
voltages, we are able to form approximately symmetric 
dots. With the relevant capacitances defined, we obtain 
the interaction energy between the two dots, using Em = 
{eyCm){{CJ:^left■C^,r^ght/Cl) - l)'' -500 /icV [22 . 
While the current structure enabled the formation of two 
nearly identical dots, our group is developing a multi (3) 
layer structure, where top-gates independently control 
the islands, a second layer of gates provides contacts to 
in-diffused source and drain, and a third layer provides 
the barrier gates. This structure allows the source and 
drain reservoirs to remain populated even for low electron 
numbers in the dots. 

Figures 3(b,c) show fine scans of bias triangles at 
VsD=^'^ and 0.5 mV respectively. Resonant tun- 
neling through the ground state and excited states of 
the double dot is clearly observed in the high-resolution 
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FIG. 3: Bias spectroscopy of a weakly coupled double dot with 
yBM=0.802 V. (a) At finite source-drain bias, the triple points 
develop into triangle pairs. Relevant capacitances of the DQD 
can be extracted from the sizes of the hexagon and triangles. 
The fractional splitting, F ^ 0.2 indicates a dominant inter- 
dot capacitive coupling, (b, c) Detailed bias spectroscopy of 
a pair of triangles at Vsd=1-0 mV and 0.5 mV. (d) Line cut 
along the red arrow in (b) shows resonant tunneling through 
excited states in the transport. 



TABLE I: Comparison of experimental values obtained from 
Fig. 3(a) and modelled FastCap capacitances. For definitions, 
see text. 
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bias-spectroscopy. With increasing Vsd^ the triangular 
conducting regions become larger with more discrete lev- 
els in the bias window and the overlap of the triangle 
pairs increases. Figure 3(d) shows a plot of Isd as a 
function of detuning energy, e ^ between levels of the 
double dot. This Isd line trace is extracted from a cut 
of the bias triangle as shown in Fig. 3(b), where the 
ground and excited state resonances are indicated by the 
labels a-d. The energy splitting of the first excited state 
b to its ground state a is ~300 /ieV. We roughly estimate 
the average energy-level spacing of a dot via Weyl's for- 
mula, AE = / gm* where A is the area of the dot. 
For a 2DEG system in Si, the effective mass of the elec- 
trons m*=0.19me and the degeneracy ^=4, taking into 
account the spin and valley degeneracies [24 . Using this 
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we calculate AE ~400 /ieV, which would be the expected 
average level spacing if all symmetries are broken. Since 
no field is applied to the dots, the spacing would be a 
factor of 2 larger or ~800 jaeV. In Fig. (d), we moni- 
tor transport through a serial configuration of two dots 
along the line cut presented in (b). In this case, we move 
the energy levels in both dots in opposite direction with 
respect to each other [2j with results in an effective re- 
duction of a factor of two in the experimentally expected 
level splitting, in good agreement with the experimental 
data. 

In conclusion, we have presented a tunable double- 
gated DQD defined in intrinsic silicon. The fabrication 
of the device is reproducible and MOS-compatible, en- 
abling scale- up or integration into more complex designs. 
Transport measurements have been performed and ex- 
tracted device capacitances were in good agreement with 
FastCap modelling. High resolution bias spectroscopy 
of the double dot presented evidence of resonant tunnel- 
ing through ground and excited states, indicating that 
the system was in the few-electron regime. To reduce 
the electron number to a single electron in each dot we 
propose the incorporation of additional plunger gates, in- 
dependently controlling each dot, together with an inte- 
grated charge detector to monitor the dot occupancies. 
Such Si-based double quantum dot structures would have 
excellent potential for the investigation of the singlet- 
triplet two-level system due the long spin-coherence times 
in silicon. 
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